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Background: While it is known that elevated baseline intracardiac repolarization lability is 
associated with the risk of fast ventricular tachycardia (FVT)/ventricular fibrillation (VF), the 
effect of its longitudinal changes on the risk of FVT/VF is unknown. 

Methods and Results: Near-field (NF) right ventricular (RV) intracardiac electrograms 
(EGMs) were recorded every 3-6 months at rest in 248 patients with structural heart 
disease [mean age 61. 2± 13.3; 185(75%) male; 162(65.3%) ischemic cardiomyopathy] 
and implanted cardioverter-defibrillator (ICD) or cardiac resynchronization therapy 
defibrillator (CRT-D) [201 (81 %) primary prevention]. Intracardiac beat-to-beat QT variability 
index (QTVInf) was measured on NF RV EGM. During the first study phase (median 18 
months), participants made on average 2.4 visits. Then remote follow-up was continued 
for an additional median period of 3 years. Average QTVInf did not change during the first 
year after ICD implantation (-0.342 ± 0.603 at baseline vs. -0.262 ± 0.552 at 6 months 
vs. -0.334 ± 0.603 at 12 months); however, it decreased thereafter (-0.510 ± 0.603 at 18 
months; P = 0.042). Adjusted population-averaged GEE model showed that the odds of 
developing FVT/VF increased by 75% for each 1 unit increase in QTVInf- (OR 1.75 [95%CI 
1.05-2.92]; P = 0.031). However, individual patient-specific QTVInf trends (increasing, 
decreasing, flat) varied from patient to patient. For a given patient, the odds of developing 
FVTA/F were not associated with increasing or decreasing QTVInf over ti me [OR 1-27; 
(95%CI 0.05-30.10); P = 0.881]. 

Conclusion: While on average the odds of FVTA/F increased with an increase in QTVInf, 
patient-specific longitudinal trends in QTVI N f did not affect the odds of FVTA/F 



Keywords: intracardiac electrograms, repolarization lability, ventricular tachyarrhythmia, longitudinal analysis, QT 
variability index 



Augmented intracardiac repolarization lability predicts ventric- 
ular arrhythmia (Haigney et al., 2004; Couderc et al., 2007; 
Tereshchenko et al, 2009), sudden cardiac death (SCD) (Piccirillo 
et al., 2007) and cardiovascular mortality (Tereshchenko et al, 
2012) in patients with structural heart disease and long QT 
syndrome (Hinterseer et al, 2008). Underlying mechanisms of 
increased repolarization lability are associated with arrhythmo- 
genic substrate [scar, fibrosis, local heterogeneities in action 
potential duration and morphology, cell-to-cell uncoupling, 
stochastic gating of ion channels (Pueyo et al., 2011), and 
increased sympathetic tone in the ventricles of the heart (Baumert 
et al, 2011)]. It is known that arrhythmogenic substrate does 
change over time. However, little is known about neither the rate 
of changes in repolarization lability over 6-12-18 month period 
nor the predictive value of these changes. 

Dynamic time-dependent association between alternating 
repolarization variability and cardiac death has been previ- 
ously reported. The Alternans Before Cardioverter Defibrillator 



(ABCD) Trial showed that baseline microvolt T-wave alter- 
nans (TWA) was significantly associated with sustained ven- 
tricular tachycardia (VT)/ ventricular fibrillation (VF) or SCD 
at 6 months, but lost association at 12 months of follow-up 
(Costantini et al., 2009). ABCD trial results brought up discus- 
sion regarding the necessity for establishing the "expiration date" 
for TWA results, and raised the question about possible utility of 
repetitive TWA assessment. While elevated repolarization lability, 
TWA and other non-invasive risk markers (including decreased 
heart rate variability) have been found to be associated with life- 
threatening VT/VF and SCD, none of them are routinely used for 
risk stratification of primary prevention implantable cardioverter 
defibrillator (ICD) (Tereshchenko and Berger, 2012). Repeated 
measurement of risk markers might strengthen their predictive 
value. However, studies of longitudinal changes in repolarization 
lability are limited. 

Earlier we showed that repolarization lability is present 
throughout the myocardium (Tereshchenko et al., 2009). We 
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demonstrated that the single baseline measurement of intracar- 
diac QT variability index (QTVInf), was associated with fast ven- 
tricular tachycardia (FVT)/VF during the subsequent 16 months 
of follow-up. However, the effect of longitudinal changes in 
intracardiac repolarization lability on the risk of FVT/VF is 
unknown. The goal of the present study was to determine if 
patients with increasing intracardiac QTVInf over time expe- 
rience a greater probability of having FVT/VF than those with 
decreasing or stable intracardiac QTVInf- 

METHODS 

The ICD-EGMs study (Tereshchenko et al, 2009) protocol was 
approved by the Johns Hopkins University and the Washington 
University Human Studies Committees and all participants gave 
written informed consent before entering the study. 

STUDY POPULATION 

The ICD-EGMs study design (NCT00916435) was previously 
described (Tereshchenko et al., 2009). Briefly, patients with struc- 
tural heart disease of either sex older than 18 years were enrolled 
in the study if they had a transvenous ICD or a cardiac resynchro- 
nization therapy defibrillator (CRT-D) device implanted for pri- 
mary or secondary prevention of SCD. In this study we included 
only participants who (1) had ICD or CRT-D implanted 1 week 
before enrollment, (2) made at least 2 consecutive office follow- 
up visits during the first study phase (2005-2007), (3) had at 
least 2 analyzable digital recordings of intracardiac EGMs at rest 
that were at least 30 days apart, (4) had EGMs recorded in sinus 
rhythm, with (5) identical type of ventricular activation [either 
ventricular-sensed (VS), or ventricular-paced (VP)] at baseline 
and at all follow-up visits. 

RECORDING OF INTRACARDIAC EGMs AND FOLLOW-UP OF STUDY 
PARTICIPANTS 

Recording of intracardiac EGMs was performed at every office 
follow-up visit (every 3-6 months) during the first 2 years of 
the study (2005-2007). Patients were then followed remotely 
via Carelink® and Latitude® during next 4 years (2007-2011). 
Near-field (NF) and far-field (FF) RV intracardiac EGMs were 
recorded at rest for 5-15 min simultaneously with surface electro- 
cardiogram (ECG) via ICD programmer as previously described 
(Tereshchenko et al, 2009). 

Duration of time periods between follow-up visits varied from 
patient to patient, due to the observational nature of the study. In 
order to standardize assessment of longitudinal changes in EGM 
parameters, we categorized follow-up periods as the following. All 
follow-up visits that occurred in a period of 1-6 months after the 
baseline (Visit 1) were considered as Visit 2. Visit 3 EGM record- 
ing was performed at any time in a period of 181-365 days after 
baseline EGM recording. Visit 4 was performed during the 1st half 
of the 2nd year of follow-up. 

INTRACARDIAC REPOLARIZATION LABILITY ANALYSIS 

Intracardiac repolarization lability was measured on NF RV 
EGM as previously described (Tereshchenko et al, 2009) by 
Berger's method (Berger et al., 1997), using customized MATLAB 
(MathWorks, Inc., Natick, MA) software. The R-wave peak 



was automatically detected on NF EGM channel. Two investi- 
gators (AG, JL) defined an intracardiac intervals template by 
selecting the beginning and the end of major NF EGM deflec- 
tion, and the end of the T wave. The algorithm then determined 
how much the repolarization segment of each beat had to be 
stretched or compressed to match the template. Appropriate 
selection of fiducial points was verified by the third investigator 
(LGT). Premature ventricular and atrial beats with one post- 
ectopic sinus beat were excluded from the analysis. Recordings 
with more than 15% of ectopic, or noise-distorted beats were 
excluded. 

ENDP0INTS 

Appropriate ICD shocks for FVT/VF served as the endpoints 
in this study. Programming of the ICD device was based on 
clinical evaluation of the attending electrophysiologist. ICD inter- 
rogation data was adjudicated by an endpoint committee com- 
posed of 3 members. ICD therapy occurring for VT or VF was 
classified as appropriate (Tereshchenko et al, 2009). FVT/VF 
was defined as VT/VF with an average cycle length (CL) < 
240 ms. After the 1st FVT/VF event follow-up was continued, 
and all subsequent sustained FVT/VF events with appropri- 
ate ICD shocks were included in the analysis as the study 
endpoints. 

STATISTICAL ANALYSIS 

The association of baseline clinical characteristics and ECG 
parameters with the type of presenting rhythm was measured by 
the x. 2 -tests and AN OVA or f-test, respectively, for categorical and 
continuous variables with normal distribution. A test of equality 
of medians or Wilcoxon rank-sum test was used in case the dis- 
tribution of parameter was not normal. A P-value of <0.05 was 
considered significant. Data analysis was performed using STATA 
12.1 (StataCorp LP, College Station, TX). 

As office visits were unequally spaced in time, QTVI vari- 
ogram was explored, and follow-up time periods were equalized 
with incremental 6 months-intervals. Possible effects of VP on 
longitudinal QTVInf changes were explored in patients with dif- 
ferent devices types (single-chamber ICD, dual-chamber ICD, 
CRT-D), as well as in patients who predominantly (at least 
99% of follow-up time and 100% of recorded EGM) had VS, 
in comparison to VP rhythm. "Spaghetti" plots were examined 
to study variations of QTVInf across time for each individ- 
ual. Smoothed plots were used to explore group response of 
QTVInf as a function of time, to study variations of QTVInf 
across different individuals. Average changes in QTVInf per 
6 months of follow-up were determined. Patterns of QTVInf 
changes over time in patients with and without FVT/VF were 
compared. Since repeated measures made on the same sub- 
ject are correlated, within-person correlations matrix of QTVInf 
across multiple visits was estimated, and the correlation struc- 
ture was described. The correlation structure for the resid- 
uals was explored after removing the mean trend effect. To 
standardize QTVInf at each visit, we subtracted observations 
in each category (each visit) by the mean for that visit and 
then divided them by the standard deviation for that visit. 
Further longitudinal regression analyses took into account the 
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visit-to-visit QTVInf correlation structure to obtain valid infer- 
ences. 

We compared the association of FVT/VF events with pre- 
ceding longitudinal QTVInf changes in an average study par- 
ticipant (in the population-averaged model), and in specific 
study subject (in the subject-specific model). Population- 
averaged marginal model accounting for correlation struc- 
ture [Generalized Estimating Equations (GEE) model] was 
developed. Multivariate GEE model was used to determine 
an association between a population-averaged longitudinal 
QTVInf changes, and a subsequent outcome. FVT/VF events 
with appropriate ICD shocks that occurred after respective 
QTVInf measurements served as an outcome. GEE model 
was adjusted by age, sex, race, left ventricular ejection frac- 
tion (LVEF), New York Heart Association (NYHA) heart failure 
class (dichotomized as NYHA = II), indication for ICD (pri- 
mary or secondary prevention of SCD), type of cardiomyopathy 
(ischemic or non-ischemic), history of revascularization proce- 
dure (either PCI or CABG), type of implanted device(single- 
chamber ICD, dual-chamber ICD, CRT-D), type of analyzed 
rhythm (VS or VP), and use of class III anti- arrhythmic medi- 
cations. 

In order to determine the patient-specific association of lon- 
gitudinal QTVInf changes with subsequent FVT/VF, we ran a 
random intercept model for FVT/VF. We accounted for corre- 
lation of the repeated QTVInf observations by including a ran- 
dom intercept for each patient and control for patient's QTVInf 



(centered). Adequate fitting of the model was checked to ensure 
that the specified quadrature has adequately approximated the 
likelihood. 

RESULTS 

STUDY POPULATION 

We analyzed the data of 248 study participants: 185 (74.6%) 
men, 201 Whites (81%), mean age 61.2 ±13.3 y. Clinical 
characteristics of study participants are shown in Table 1. 
We categorized study participants based on their presenting 
rhythm during follow-up office visits. VS participants were 
paced <1% during the study duration (Table 1) and were not 
paced from ventricle(s) during study EGM recordings. Patients 
who presented with a sinus rhythm, and were 100% VP dur- 
ing the study EGM recordings, comprised VP group. Figure 1 
shows representative examples of longitudinal changes in RR' 
and QT variability in patients with and without ventricular 
arrhythmia. 

As expected, due to the differences in HF severity, CRT-D 
patients had significantly higher mean heart rate as compared to 
single-and dual-chamber ICD patients (Table 2). However, QT 
interval on NF RV EGM (QTnf) was longer in patients with dual- 
chamber ICDs as compared to those with single chamber ICDs 
and CRT-Ds. There were no significant differences in heart rate 
variance, QTnf variance, and QTVInf amongst patients with dif- 
ferent device types, nor in patients with presenting VS vs. those 
with VP sinus rhythm (Table 2). 



Table 1 | Comparison of demographic and clinical characteristics of patients with single-, dual-chamber ICD, and CRT-D, and in patients with 
sinus ventricular-sensed and ventricular-paced rhythm. 





Single chamber 
ICD (n = 139) 


Dual chamber 
ICD (n = 97) 


CRT-D [n = 12) 


ANOVA P 


V-sensed 
(n= 195) 


V-paced 
(n = 53) 


P-value 


Age(SD), y 


58.2 (13.0) 


64.5 (13.3) 


69.8 (7.5) 


0.0001 


60.3 (13.4) 


64.5 (12.6) 


0.036 


Males, n(%) 


101 (72.7) 


74 (76.3) 


10 (83.3) 


0.636 


145 (74.4) 


40 (75.5) 


0.869 


African American, n(%) 


35 (25.2) 


10 (10.3) 


2 (16.7) 


0.016 


40 (20.5) 


7 (13.2) 


0.229 


CHF NYHA class = II, n(%) 


92 (66.2) 


52 (53.6) 


12 (100.0) 


0.003 


118 (60.5) 


38 (71.7) 


0.135 


Ischemic CM, n(%) 


85 (61.2) 


69 (28.9) 


4 (33.3) 


0.283 


70 (35.9) 


16 (30.2) 


0.439 


Primary prevention of SCD, n(%) 


28 (20.1) 


19 (19.6) 


0 


0.228 


157 (80.5) 


44 (83.0) 


0.680 


Presenting VP rhythm, n(%) 


2 (1 .4) 


40 (41.2) 


11 (91.7) 


<0.0001 








Atrial pacing %, median(IQR) 




11.5 (0-76) 


0 (0-0) 


<0.0001 


16 (0-76) 


1.5 (0-74) 


0.328 


Ventricular pacing %, median(IQR) 


0 (0-0) 


4.5 (0-29.5) 


100 (98-100) 


<0.0001 


0 (0-1) 


16 (1-91) 


<0.0001 


LVEF(SD), (%) 


31.9 (12.2) 


32.4 (11.3) 


25.0 (12.4) 


0.127 


32.4 (12.1) 


29.3 (11.2) 


0.075 


Diabetes mellitus Hx, n(%) 


50 (36.0) 


23 (24.0) 


3 (25.0) 


0.132 


65 (33.5) 


11 (20.8) 


0.075 


Hypertension Hx, n(%) 


108 (77.7) 


67 (69.8) 


10 (83.3) 


0.306 


150 (77.3) 


35 (66.0) 


0.093 


CABG or PTCA Hx, n(%) 


65 (46.8) 


51 (52.6) 


7 (58.3) 


0.561 


94 (48.2) 


29 (54.7) 


0.400 


History of AF, n(%) 


32 (23.0) 


43 (44.3) 


7 (58.3) 


<0.0001 


55 (28.2) 


27 (50.9) 


0.002 


Beta blockers, n(%) 


126 (90.7) 


78 (80.4) 


8 (72.7) 


0.038 


171 (87.7) 


41 (78.9) 


0.104 


ACE-lnhibitors or ARBs, n(%) 


111 (79.9) 


78 (80.4) 


11 (91.7) 


0.609 


155 (79.5) 


45 (84.9) 


0.376 


Digitalis, n(%) 


45 (32.4) 


39 (40.2) 


5 (45.5) 


0.375 


69 (35.4) 


20 (38.5) 


0.681 


Statin, n(%) 


97 (69.8) 


68 (70.1) 


8 (72.7) 


0.979 


137 (70.3) 


36 (69.2) 


0.886 


Nitrates, n(%) 


36 (26.1) 


20 (20.6) 


5 (45.5) 


0.170 


50 (25.8) 


11 (21.2) 


0.493 


Aldosterone antagonists, n(%) 


53 (38.1) 


37 (38.1) 


5 (45.4) 


0.888 


76 (39.0) 


19 (36.5) 


0.748 


Antidepressants, n(%) 


30 (21.6) 


33 (34.0) 


3 (27.3) 


0.105 


47 (24.1) 


19 (36.5) 


0.072 


Class III antiarrhythmics, n(%) 


24 (17.3) 


33 (34.0) 


5 (45.5) 


0.004 


42 (21.5) 


20 (38.5) 


0.012 
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FIGURE 1 | Representative example of longitudinal changes in beat-to-beat intracardiac QT and RR' variability in patient with (A) and without (B) 
sustained FVT/VF and appropriate ICD shocks. 



Table 2 | Comparison of all-visits averaged ECG parameters in patients with single-, dual-chamber ICD, and CRT-D, and in patients with sinus 
ventricular-sensed and ventricular-paced rhythm. 





Single chamber 


Dual chamber 


CRT-D (n = 24) 


ANOVA P 


V-sensed 


V-paced 


P-value 




ICD (n = 332) 


ICD [n = 230) 






[n = 464) 


(n = 122) 




Heart rate(SD), bpm 


74.2 (14.5) 


71.1 (13.3) 


78 (17.5) 


0.009 


72.8 (14.3) 


74.5 (14.2) 


0.231 


QT interval(SD), ms 


469 (116) 


496 (110) 


478 (180) 


0.030 


479 (116) 


487 (123) 


0.477 


Heart rate variance, 


40.2 (7.9-91.5) 


47.1 (12.8-99.4) 


45.7 (14.0-117.7) 


0.549 


41.0 (9.1-93.2) 


48.6 (13.0-106.7) 


0.209 


median(IQR), ms 2 
















QT variance, 


748(122-1693) 


729 (167-1704) 


556 (269-1190) 


0.699 


726 (130-1699) 


715 (196-1678) 


0.617 


median(IQR), ms 2 
















QTVI(SD) 


-0.307 (0.548) 


-0.402 (0.648) 


-0.293 (0.641) 


0.163 


-0.344 (0.570) 


-0.345 (0.679) 


0.983 



VENTRICULAR TACHYARRHYTHMIA DURING FOLLOW-UP observed in 20 (77%) out of 26 patients, frequently in clusters. 

During the first study phase participants made on average 2.4 vis- Baseline clinical characteristics in patients with or without 

its. FVT/VF events with appropriate ICD shocks were diagnosed FVT/VF during follow-up did not differ (Table 3). Importantly, 

in 26 (10.5%) patients. Multiple FVT/VF during follow-up were QTVInf was the only ECG parameter that differentiated patients 
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Table 3 | Comparison of baseline clinical and demographic 
characteristics in patients with and without FVT/VF during follow-up. 



Age(SD), y 

Males, n(%) 

African American, n(%) 

CHF NYHA class = II, n(%) 

Ischemic CM with Ml history, 

n(%) 

Primary prevention of SCD, 
n(%) 

Single-chamber ICD, n{%) 
Dual-chamber ICD, n(%) 
Bi-Ventricular ICD, n(%) 
Atrial pacing 1 %, 
median(IQR) 
Ventricular pacing 2 %, 
median(IQR) 
LVEF(SD), (%) 
Diabetes mellitus, n(%) 
Hypertension, n(%) 
CABG or PTCA, n(%) 
History of AF, n(%) 
Beta blockers, n{%) 
ACE-lnhibitors or ARBs, n(%) 
Digitalis. n(%) 
Statin, n(%) 
Nitrates, n(%) 
Aldosterone antagonists, 
n(%) 

Antidepressants, n(%) 
Class III antiarrhythmic 
medication, n(%) 



Percentage of atrial pacing during study phase ? (counters data) in patients 
with dual-chamber ICD or CRT-D; 2 percentage of ventricular pacing, accordingly, 
in patients with ICD only (CRT-D excluded). 

with and without arrhythmia. QTVTnf was significantly higher in 
patients with FVT/VF than in patients without FVT/VF (Table 4). 

THE MEAN TREND IN LONGITUDINAL CHANGES OF INTRACARDIAC 
REPOLARIZATION LABILITY 

Average QTVInf in study participants with single-chamber ICD 
did not change during the study period. However, average 
QTVTnf in patients with dual-chamber ICDs increased dur- 
ing the first year after ICD implantation, but slightly decreased 
thereafter (Figure 2). At the same time, no statistically signifi- 
cant differences were observed in all VS, as compared to all VP 
patients (Figure 3). Averaged trend in QTVTnf was different in 
patients with FVT/VF, as compared to patients without FVT/VF 
(Figure 4). 

SUBJECT-SPECIFIC CHANGES IN REPOLARIZATION LABILITY 
OVER TIME 

Individual patients demonstrated drastically different trends: in 
some patients intracardiac QTVInf increased while in others 



No FVT/VF 


FVT/VF 


P-value 


[n = 222) 


(n = 26) 




61.9 (12.8) 


55.3 (16.6) 


0.058 


165 (74.3) 


20 (76.9) 


0.773 


43 (19.4) 


4 (15.4) 


0.624 


138 (62.3) 


18(69.2) 


0.480 


145 (65.3) 


17 (65.4) 


0.994 


180 (81.1) 


21 (80.8) 


0.969 


124 (55.9) 


15 (57.7) 


0.963 


87 (39.2) 


10 (38.5) 


0.963 


11 (5.0) 


1 (3.9) 


0.963 


9.5 (0-71.5) 


32 (1-92) 


0.053 


0 (0-3) 


0 (0-0) 


0.152 


31.7 (11.8) 


32.5 (13.5) 


0.777 


68 (30.8) 


8 (30.8) 


1.00 


165 (74.7) 


20 (76.9) 


0.901 


111 (50.0) 


12 (46.2) 


0.711 


74 (33.3) 


8 (30.8) 


0.793 


192 (86.9) 


20 (76.9) 


0.169 


181 (81.5) 


19 (73.1) 


0.302 


76 (34.4) 


13 (50.0) 


0.117 


156 (70.6) 


17 (65.4) 


0.584 


55 (25.0) 


6 (23.1) 


0.830 


86 (38.9) 


9 (34.6) 


0.670 


58 (26.2) 


8 (30.8) 


0.622 


55 (24.9) 


7 (26.9) 


0.821 



QTVInf decreased, stayed flat, or fluctuated over time. There were 
no differences in the behavior of the subject-specific longitudinal 
relationships in patients with different device type (Figure 5), or 
any between groups of VS and VP patients (Figure 6). Moreover, 
no differences in subject-specific QTVInf longitudinal trends 
were observed in patients with FVT/VF, as compared to patients 
without FVT/VF (Figure 7). Multiple cross-overs of the lines 
indicated that the relative order of patients, ordered by their 
baseline QTVInf, changed over time. The study population- 
averaged QTVInf longitudinal relationship (Figures 2-4) was 
not consistent with the subject-specific longitudinal relationship 
(Figures 5-7). 

We observed weak positive correlation between QTVInf at 
the 1st and the 4th visits, and at the 1st and the 2nd visits 
in patients with single-chamber ICD (r = 0.312; P = 0.047 and 
r = 0.269; P = 0.024, respectively) and VS patients (r = 0.297; 
P = 0.026 and r = 0.384; P < 0.0001, respectively), and weak 
negative correlation in QTVInf at the 3rd and 4th visit (r = 
-0.604; P = 0.003 in single-chamber ICD group; r = -0.426; 
P = 0.021 in VS patients). Thus, in a given patient with single- 
chamber ICD who did not experience ventricular pacing, QTVInf 
1 week after ICD implantation more likely positively corre- 
lated with QTVInf 1.5 years after ICD implantation. If such 
a patient experienced elevation of QTVInf during the 1st year 
post-ICD implantation, then during the next 6 months QTVInf 
was more likely decreasing (negative correlation in QTVInf 
between the 3rd and the 4th visits). No significant correla- 
tions between QTVInf observations at different visits in VP 
patients, both with dual-chamber ICDs and CRT-D devices were 
found. 

ASSOCIATION BETWEEN LONGITUDINAL CHANGES IN 
REPOLARIZATION LABILITY AND SUBSEQUENT FVT/VF 

The mean QTVInf trend in patients without FVT/VF 
demonstrated slight, but significant decrease 1.5 years after 
device implantation (Figure 4A), whereas no changes in 
mean QTVInf were observed in patients with FVT/VF 
(Figure 4B). Patterns of the subject-specific relationships 
between QTVInf and time looked alike in patients with and 
without FVT/VF (Figures 7A,B). QTVInf correlations structure 
in patients without FVT/VF revealed weak positive correla- 
tion between the 1st and the 2nd visit (r = 0.257; P = 0.004) 
and negative correlation of approximately the same strength 
between the 3th and the 4th visit (r = -0.339; P = 0.040). 
However, QTVInf measured at different visits in a given 
patient who experienced FVT/VF during follow-up, did not 
correlate. 

In order to study patient-specific dynamic changes in intracar- 
diac RL before FVT/VF, we plotted QTVInf before each FVT/VF 
event. We used the actual time from then EGM recording to the 
FVT/VF event as a continuous variable (Figure 8). Noticeably, 
QTVInf distribution in patients without FVT/VF (Figure 8A) 
looked similar to that in patients before FVT/VF (Figure 8B). 
Trends of increasing, decreasing, and flat over time QTVInf were 
observed before FVT/VF events (Figure 8C). Of note, consis- 
tent pattern of increasing over time QTVInf before all FVT/VF 
events was observed in some (but not all) individual patients 
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Table 4 | Comparison of all-visits averaged ECG parameters in patients with and without FVT/VF during follow-up. 



No FVT/VF 
(n = 525 visits) 



FVT/VF 

(n = 61 visits) 



P- value 



Heart rate(SD), bpm 
QT interval(SD), s 

Heart rate variance, median(IQR), ms 2 
QT variance, median(IQR), ms 2 
QTV (SD) 



73.1 (14.0) 
0.479 (0.116) 

45.1 (10.4-99.4) 
697.0 (147.0-1642.3) 
-0.362 (0.601) 



73.4(16.3) 
0.498 (0.129) 
31.9 (12.8-70.4) 
1008 (175.0-2022.4) 
-0.186 (0.512) 



0.899 
0.278 
0.143 
0.238 
0.015 



A , Single-chamber ICD 



- 038 (0.52) 

t 



i 7 



1 
1 



ANOVAP = 0,120 



j Dual-chamber ICD 


-0.46(0.66) -0-29(0-58) -0.26(0.68) 




\ -0.63(0.61) 

¥ \ 
| \ 


\ 


f \ 


> 




* 

ANOVAP = 0.031 



CRT-D 

-0.3810.93) 



FIGURE 2 | Comparison of population-averaged longitudinal trends in QTVInf in patients with (A) single-chamber ICD, (B) dual-chamber ICD, (C) 
CRT-D. Scatterplot of QTVI NF at each standardized study visit is shown. Lowess smoother curve shows mean trend in QTVI NF over time. 



Ventricular - sensed b 



Ventricular - paced 



- 0.32 (0.59) 



- 0.56 (0.68) 



ANOVAP = 0.163 



ANOVAP = 0.1 



FIGURE 3 | Comparison of population-averaged longitudinal trends in QTVInf in patients with (A) presenting ventricular-sensed (VS), and (B) presenting 
ventricular-paced (VP) sinus rhythm. Scatterplot of QTVInf at each standardized study visit is shown. Lowess smoother curve shows mean trend in QTVInf over time. 



with multiple FVT/VF events. In contrast, patients with a sin- 
gle FVT/VF event tended to demonstrate rather decreasing over 
time QTVInf. However, small subgroups size did not allow us to 
quantify observed differences. 

GEE POPULATION-AVERAGED MODEL 

In multivariate GEE analysis with independent correlation struc- 
ture, increasing over a time course of several months QTVInf 
was significantly associated with FVT/VF event [OR 1.75 (95%CI 
1.05-2.92); P = 0.031]. Thus, on average, for patients in our 
study, the odds of developing FVT/VF increased by 75% for each 
1 unit increase in QTVInf- 



SUBJECT-SPECIFIC RANDOM INTERCEPT MODEL 

Association between patient-specific dynamic changes in QTVInf 
and FVT/VF was studied in the random-effects logistic regression. 
In order to determine if a specific patient with increasing QTVInf 
over time experiences a greater probability of having FVT/VF, 
as compared to patient with decreasing or stable QTVInf, 
we ran random-effects logistic regression analysis. For a given 
patient, the odds of developing FVT/VF were not associated with 
increasing or decreasing QTVInf [OR 1.27; (95%CI 0.05-30.10); 
P = 0.881] over time. We used 194 integration points in this 
model for assurance that the likelihood is appropriately approxi- 
mated, which was confirmed. In addition, we explored potential 
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-0.37(0.62) -0.27(0.551 



No FVT/VF 

- 035 (0.59) 



•1 

f-H 



ANOVA P=0.051 



FVT/VF 



ANOVA P = 0.941 



FIGURE 4 | Comparison of population-averaged longitudinal trends in QTVInf in patients (A) without FVT/VF, and (B) with sustained FVT/VF and 
appropriate ICD shocks. Scatterplot of QTVI NF at each standardized study visit is shown. Lowess smoother curve shows mean trend in QTVI NF over time. 




FIGURE 5 | "Spaghetti" plots of individual patient-specific longitudinal relationships between QTVI NF and time for each subject with (A) 
single-chamber ICD, (B) dual-chamber ICD, (C) CRT-D. 



Ventricular - sensed 



Ventricular - paced 




FIGURE 6 | "Spaghetti" plots of individual patient-specific longitudinal 
relationships between QTVI NF and time for each subject with (A) 
presenting ventricular-sensed (VS), and (B) presenting 
ventricular-paced (VP) sinus rhythm. 




FIGURE 7 | "Spaghetti" plots of individual patient-specific longitudinal 
relationships between QTVInf and time for each subject (A) without 
FVT/VF, and (B) with sustained FVT/VF and appropriate ICD shocks. 



interactions and adjusted patient-specific models for the effect 
of VP in preceding follow-up interval, for percentage of beats, 
included in QTVI analysis by automated software, and for the 
differences in follow-up time interval. Subject-specific longitu- 
dinal QTVInf trend was associated neither with study outcome 
(FVT/VF), nor with any other measured in this study clinical, 
demographic, or ECG parameter. 

DISCUSSION 

Our study revealed differences in population-averaged and 
patient-specific intracardiac QTVInf trends over 1.5 years 
after device implantation. In patients without FVT/VF events 
population-averaged QTVInf decreased after the 1st year since 
ICD implantation, whereas in patients with FVT/VF events 



averaged QTVInf did not change. Group-averaged QTVInf 
was independently associated with the odds of life-threatening 
FVT/VF. On average, the odds of FVT/VF increased by 75% for 
each 1 unit increase in QTVInf. In contrast, individual subject- 
specific QTVInf trends in many patients differed from the mean 
trend for the whole study population. Individual subject-specific 
trends in intracardiac QTVInf were not associated with subse- 
quent FVT/VF. 

LONGITUDINAL STUDIES OF REPOLARIZATION LABILITY: MEAN TREND 
VS. SUBJECT-SPECIFIC CHANGES OVER TIME 

Population-averaged trends in repolarization lability and TWA 
changes immediately before VT/VF have been previously studied. 
The vast majority of investigators explored short-term (24-48 h) 
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FIGURE 8 | (A) Variogram of QTVI NF in patients without FVTA'F shows actual 
time since EGM recording till censored. (B) Variogram of QTVInf in patients 
with FVTA'F shows actual time since EGM recording till sustained FVTA/F 



with ICD shock. (C) Line "spaghetti" plots of the longitudinal relationships 
between QTVInf and actual time of FVTA/F for each patient with FVTA/F 
event. 



dynamic ECG changes. Shusterman et al. (2006) described an 
upsurge of alternating and non-alternating RL lOmin before the 
onset of VT/VF in HF patients. Analysis of 24-h Holter ECG 
revealed pronounced diurnal variations in QTVI (higher QT 
variability in the morning and during the day, and lower QT vari- 
ability during the night) in HF patients (Dobson et al., 2009). 
At the same time, no significant changes in the mean QTVI 
trend was detected during 12 h of ECG monitoring before the 
onset of spontaneous VT/VF in the acute intensive cardiac care 
unit patients (Sachdev et al., 2010). In this study population- 
averaged model reflected absolute difference in QTVInf 12-18 
month after ICD implantation in patients with and without sub- 
sequent FVT/VF (Table 4; Figure 4), and re-confirmed our previ- 
ous finding of association between elevated QTVInf and VT/VF 
(Tereshchenko et al, 2009). 

Very few investigators studied subject-specific trends in 
repolarization lability. Swerdlow and co-investigators (2011) 
used subject-specific longitudinal analysis approach and showed 
increased intracardiac TWA/variability immediately before 
VT/VF onset in ICD patients. To the best of our knowl- 
edge, our study is the first that compared patient-specific 
and population-averaged longitudinal trends in intracardiac 
QTVInf- Mechanistically, each individual's longitudinal out- 
come is governed by subject-specific disease dynamics over 
time. Hence, each subject's repolarization lability trajectory 
neither necessarily progresses in accord with a rigid population 
mean curve nor varies with bounded variation. Evidently, 
some "latent" factors (not measured in this study) impacted 
individual subject-specific trends in QTVInf- In patients with 
structural heart disease and systolic HF multiple parame- 
ters might change over time: coronary perfusion, myocardial 
contractility and compliance, use of medications, affecting 
repolarization (Tereshchenko et al., 2009), level of physical 
activity and autonomic balance, kidney function, and many 
others. Any of these unmeasured in this study factors might 
play a major role in the observed subject-specific QTVInf 
changes. 

Recently we showed that both high and low repolarization 
lability is associated with the risk of SCD in general pop- 
ulation (Tereshchenko et al, 2012). In our previous analysis 
of ICD-EGMs study we observed that QTVInf was decreased 



after premature ventricular contraction. Paradoxically, decreased, 
but not increased QTVInf was associated with VT/VF (Das 
et al, 2012). While arrhythmogenesis, associated with elevated 
repolarization lability is well understood (Tereshchenko and 
Berger, 2011), observation of association between decreasing 
QTVInf and subsequent FVT/VF is novel and prompts fur- 
ther investigations. As previously shown, static QT/RR rela- 
tionships (Batchvarov et al, 2002), as well as the dynamic 
pattern of QT/RR hysteresis (Malik et al, 2008) is highly 
patient-specific. It was even suggested that the individual 
QT-RR relationship has unique "finger-print-like" properties 
(Malik et al., 2008). This fact might contribute to the high 
degree inter-subject variability of the dynamic QTVInf pat- 
terns, observed in this study. Especially intriguing was obser- 
vation of consistent increasing QTVInf trend before each 
FVT/VF event in some (but not all) patients with multiple 
FTV/VF events. Further mechanistic studies are needed to explore 
possible mechanisms of arrhythmiogenesis, associated with 
decreasing QTVInf> in order to explain subject-specific trends 
fully. 

EFFECT OF VENTRICULAR PACING ON REPOLARIZATION LABILITY 

In this study we did not find significant differences in intracardiac 
QTVInf between VS and VP patients in sinus rhythm. This obser- 
vation is in concordance with our earlier study (Tereshchenko 
et al, 2011) of intracardiac QTVInf in CRT patients. Of note, 
we longitudinally studied only patients with consistently identical 
type of the presenting rhythm across multiple visits. 

Since deleterious effect of RV pacing is well known (Wilkoff 
et al., 2002), we focused our analysis on evaluation of the possi- 
ble longitudinal effects of amount of RV pacing on QTVInf and 
study outcomes. We did not find significant differences in per- 
centage of VP in patients with and without FVT/VF. Furthermore, 
the percentage of VP did not influence QTVInf and was not 
responsible for fluctuations of QTVInf from visit-to-visit. 

CLINICAL PERSPECTIVE 

In clinical practice and clinical research multiple physiolog- 
ical parameters are measured repetitively and are studied 
longitudinally. In the vast majority of longitudinal studies 
population-averaged analysis is the only analysis applied. Very few 
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longitudinal studies in electrophysiology report results of subject- 
specific longitudinal analysis. However, the physician manages the 
specific individual patient, rather than the "population-averaged" 
patient, and understanding of the patient-specific trends is 
extremely important for clinical decisions in the era of individ- 
ualized medicine. In this study, population-averaged and subject- 
specific longitudinal trend in repolarization lability demonstrated 
different degree of association with subsequent FVT/VF. Further 
prospective longitudinal studies of repolarization lability should 
be conducted in order to understand behavior and predictors of 
subject-specific longitudinal trends in repolarization lability, and 
to determine its association with clinically-meaningful outcomes. 

LIMITATIONS 

Several limitations have to be acknowledged. First of all, the EGM 
recordings were obtained not exactly at the same time of the day 
during all follow-up visits. However, previous circadian study of 
QT variability did not find significant differences in QT variabil- 
ity during the day (Dobson et al., 2009). All EGMs were recorded 



during the morning and middle-day hours (8 am-3 pm), which 
minimized possible circadian effect. 

Time intervals between follow up visits varied due to the 
observational nature of the study. In order to overcome this lim- 
itation, we standardized time intervals for longitudinal analysis. 
In addition, we employed specific analytical approach, which 
incorporated actual time from EGM recording to outcome. 
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